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SUMMARY 

A FORTRAN IV computer code for determining temperature histories in a propellant 
tank due to wall and internal heating was written. The analysis, on which the code is 
based, is presented. The assumptions and approximations made are also discussed. 

Also included are instructions for use of the code, a sample problem, a flow chart and a 
complete listing of the code. 


INTRODUCTION 

In designing space vehicles using cryogenic liquid propellants, it is necessary to know 
how wall and internal heating affect the temperature distribution within the propellant. 
Many different analytical models have been developed to predict these effects on temper- 
ature profiles. Reference 1 is a review article which discusses the different approaches 
to this problem. An analytical model developed at Lewis Research Center is discussed 
in detail in reference 2. Computer calculations based on that analysis, and presented in 
references 2 to 4 agreed well with experimental data. Recently, renewed interest in the 
problem has led to modifications and expansion of the original code based on the model of 
reference 2. The process of developing the computer code from the theory of reference 2 
is presented in this report. Because a thorough treatment of the theory exists in refer- 
ence 2, only the detail necessary to show the development of the computer code is pre- 
sented herein. Also presented are instructions for use of the resulting FORTRAN IV 
code, a sample problem, a flow chart, and a complete listing of the code. 



THEORY 


The theory presented in reference 2 was developed for determining the temperature 
profiles in an outflowing subcooled fluid subjected to both wall and internal heating. 1 
The assumed analytical flow model was based on results from small-scale experiments 
performed at Lewis. These experiments are described in reference 5. The results 
showed that, when a subcooled fluid is subjected to both nonuniform internal heating and 
wall heating, two distinct temperature regions are developed. In the lower region the 
fluid is thoroughly mixed and maintains a uniform temperature profile. In the upper 
region or stratified layer a temperature gradient is formed from the accumulation of 
warm fluid from the boundary layer along the tank walls. This is illustrated in figure 1(a) 
for a typical propellant tank. It was further observed that the temperature profiles in the 
stratified layer exhibited similarity. As used herein similarity is the property that two 
temperature profiles 0(X, t) at different times t differ only by a scale factor in X and 
9. The profile 9(X, t) is the temperature difference between T(X,t) and the initial tem- 
perature Tj, and X is the axial position within the tank. (All symbols are defined in 
appendix A. ) The analytical flow model thus assumed that the temperature profile consists 


Initial tem- Surface tem- 



(a) Heating in typical propellant tank. (b) Temperature profile for typical 

propellant tank. 

Figure 1. - Schematic diagram ol flow model. 


*The theory was developed primarily for internal heating caused by nuclear radi- 
ation. However, the theory and resulting computer code are applicable for any form of 
internal or bulk heating, for example, propellant heated by conduction from the tank 
bottom or other similar heat sources that contribute to the temperature profile in a 
manner indistinguishable from nuclear radiation heating. 
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of two parts: a lower region of constant temperature due to the internal heating and an 
upper region or stratified layer that has a temperature profile which exhibits the property 
of similarity and which is due primarily to the heat transfer from the walls. To simplify 
the analysis the following additional assumptions were made: 

(1) The fluid is subjected to a constant ullage pressure with no heat or mass transfer 
across the gas -liquid surface. 

(2) The liquid surface is at saturation temperature corresponding to the constant 
ullage pressure. 

(3) The fluid flows out of the system at a constant rate. 

(4) The heat input to the subcooled fluid originates from two sources : heat transfer 
from the tank walls and heat absorbed internally such as heat from nuclear radiation. 

(5) The heat input does not vary in the radial or circumferential directions. 

(6) The resulting temperature profiles do not vary in the radial or circumferential 
directions, except for the very thin boundary layer along the tank walls. 

(7) The physical properties of the fluid do not vary appreciably over the temperature 
range involved. These assumptions lead to the following expression for the temperatures 
in the stratified layer 


e(x, t) = e s jf(t) [l - <Mx,t)] + ^(x,t)j 


(i) 


where 



( 2 ) 


and 


<MX,t) 


e(x,t) - e b (t) 


(3) 


where 0 b (t) is the temperature difference (T b (X, t) - T.) in the uniform or bulk temper- 
ature region, 0 is the saturation temperature difference (T„ - Tj corresponding to the 
ullage pressure, and i//(X, t) is a similarity parameter and may be thought of as being 
the contribution to 9(X, t) due to wall heating. The term f (t) [l - ip(X, t)] is the necessary 
contribution of internal heating in the stratified layer to preserve the property of simi- 
larity. A typical profile illustrates the flow model in figure 1(b). Because i//(X, t) 
vanishes below the stratified layer, equation (1) can be used as the temperature profile 
in the entire fluid. The parameter \p(X, t) is assumed to have the following form 
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/X-X 0 (t)\ n 


X > X^t) 


i//(x,t) = 0 X<XQ(t) 

where XqW is the axial position of the bottom of the stratified layer, 6(t) is the thickness 
of the stratified layer, and n is a parameter which will be determined later. 

The term Xg(t) can also be written as X s (t) - 6(t) where X g (t) is the position of the 
liquid surface. The period of growth of the stratified layer, called the initial period, is 
then defined by 0 < Xq ^ X g . A later period is then defined by 0 < X_ ^ X g - where X g ^ 
is the location of the liquid surface at XQ(t) = 0. In the later period equation (1) is re- 
placed by 


where equations (1) and (4) are matched at tg corresponding to Xg(t) ,0 by the condi- 
tion that f(tg) = F(t 0 ). 

In order to solve for the unknowns, f(t), F(t), X_ , , 6(t), and n, an equation for the 

S y C 

energy balance for the system and an equation for the energy balance between the boundary 
layer and the wall heating are used. 

For convenience, the variable X 0 is introduced as the independent variable through 
the transformation. 



where L is the initial liquid level, Wp is the flow rate, p is the density of the propel- 
lant, and A(t) is the cross-sectional area of the propellant tank. 

The equations which result from the theory and their numerical solution are discussed 
in the following section. 

NUMERICAL SOLUTION OF EQUATIONS 

Of the unknowns, f(t), F(t), X g l , 6(t), and n, the equation for n is the most 
straightforward and is given by 



0(X,t) = flJF(t) 1 - [1 


X c (t) - X\ n 
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( 6 ) 


n = 4 — - 1 
e w 


where 0 W is the temperature rise across the boundary layer. 

Equation (6) was developed from the vertical-flat plate turbulent free -convection 
boundary -layer theory. 

The solutions for X^. ^ and 6(X g ) both come from the same equation which is 

j2 


r 


d 

dX_ 


6(X S )A(X S ) 


5 2 (XJ -1 
s dX. 


• _ A <V 


6 S (X ) 


[ A < x t 




(n + 1) 


(n + l)(n + 2) 


(n + l)(n+ 2)(n + 3) 


.1 


V 


A (X S ) j 
C pVs J 

•'o 


s 


q * (x) dl dx 


(7) 


where C p is the specific heat at constant pressure of the fluid, q(X) is the wall heat 
-t w 

flux, and cr is the surface area of the tank. Let 


<W = 


■/ 


q » (x) s dx 


As will be seen later Q^Xg) is essential input information for the computer code. 

To solve for X 7 , integrate both sides of equation (7) with respect to X between 
the limits L and X 7 . Noting that 6(L) = 0 leads to the following equation: 


a < X s,P A < x s, ; > 

(«+ 1 ) 


a2 < x s, ; > jj f- [ a(x m>] [ a(x m>] 

s, l , Sjl 


(n + l)(n + 2) 



(n + l)(n + 2)(n + 3) 
•X„ 


v s ,1 


MX) 

<V X s> 


dX„ 


C P W P e s 


( 8 ) 


Noting that 6(X 7 ) = X 7 gives the following equation 

S y L S J i. 
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(n + 1) 


X 2 d I 

8 > l dX J 

S, t 


(n + l)(n + 2) 


’ dX 2 L 8,1 
s, Z 

(n + l)(n + 2)(n + 3) 


L 


A «s> „ (x)dx _ 
c Ih e W X s “ + 
'■'P w P t s 


/ 


A (x s ) Q w ( x s ) 

C pVs 


dX D (9) 


Equation (9) can now be put in a form suitable for solution in the computer code by letting 
the left side be Z^(X g j) and the right side be Z R (X g ^), then 

Z L< X s,;>-V X s,!> = 0 (10) 

This equation is solved in the computer code by iteration, that is, by letting 

Z L< X s,;>- Z R< X s,!>' R < X s,i> < U) 

and continually readjusting X_ , until R(X_ , ) becomes sufficiently close to zero. The 

S y L S y L 

integration involved for determining Z R (X g p is performed within the code by using a 
straightforward Simpson's rule integration. 

Equation (7) is also used, as mentioned previously, to calculate S(X g ). Substituting 
on the right and taking the derivative on the left gives, 


6- (X ) 28(X s> A '( X s> , 352 (* s )A"(X 3 ) 

S Un + 1) (n + l)(n + 2) (n + l)(n + 2)(n + 3) 

p>(X g )A’(X g ) 5 2 (X g )A"(X g ) 6 3 (X s )A"*(X s ) 

I (n+ 1) (n + l)(n +2) (n + l)(n + 2)(n + 


A(X ) 

QuW (12) 

c pW p e s ^ 


Solving for 6’(X g ) from equation (12) and using a Taylor series expansion for 6(X g - AX g ) 
about X g , neglecting all but the first two terms, gives the expression for 6(X g - AX g ) 
used in the computer code 

6(X g - AX g ) = 6(X g ) - 6’(X g )AX g (13) 
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The remaining unknowns are calculated by solving the following equations for f(X ) 

" s 


and F(X g ) 


df(X ) A(X ) 

a(X ) ®- - b(X )f(X ) = -L-QJk ) 

s dX s s C p W p 0 n s 


(14) 


and 


dF(Xj r -i 

" a(Xs) ~ir ~ ' r(Xs) i 1 ‘ p( Vl = 


A < X s> 

C P W P*s 


%(X S ) + Q n (X g ) (15) 


'] 


where 


<W - 


b(X s ) = - 


r 


q n(X)A(X)dX 


A < X s> 

C P^P 0 s 


Qw< X s> 


(16) 

(17) 



and q n (X g ) is the nuclear or internal heating rate per unit volume. 

The second integral indicated in equation (18) for a(X g ) is approximated by the 
following equation: 
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b(X)dX = 


M 

i=0 


f b(L - iAX s ) + b[L - (i + 1)AX J 


►AX_ 


( 21 ) 


where 


M = - 1 

AX s 

Both quantities f(X g - AXg) and F(X g - AX g ) may be expanded in a Taylor scries 

about X g . Neglecting all but the first two terms of each series and solving for df/dX g 

and dF/dX. from equations (14) and (15), respectively, give the following expressions 
s 

which are used in the code : 


f(X g - AX g ) = f(X g ) + 


'W A ( X s> 

a( x g )CpW p 


bag 

a( x s ) 


f ( x s ) 


AX„ 


( 22 ) 


F(X g - AX g ) = F(X g ) 



( X s> + 


Q-(X g )C p Wp0 


A(x s ) + r(x g ) 


s 



(23) 


USING CODE 

The quantities Q w (X g ) and Q n (X g ) that appear in the previous equations are essen- 
tial input to the computer code. The values needed as a function of X g are obtained 
from two subroutines WALL and NUC which are written by the code user for the heating 
configuration of concern. 

The computer code was written specifically for a propellant tank with the geometry 
described in table I. The code may also be used for a cylindrical tank with a spherical 
bottom by letting Rq = Rj = d = Xj = Xg = Xg. However, the code may be rewritten for 
propellant tanks with other geometries. To do this, new solutions for variables which 
depend on A(X g ) must be developed. This is a straightforward process especially in 
situations where there are only two or three different geometric regions. 

The following lists contain the essential input information and a description of the 
output for the existing code. 
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Region a 

Equation of tank 
profile, R(X) 

Cross-sectional area, A(X) 

Geometry 

OsX<Xj 

[ R o- ^o-x) 2 ] 172 

jt [-X 2 + 2Rqx] 

Sphere 


X+Rq^- 1) 

jrfx 2 + 2R 0 (v^"- 1)X+ Rq(/T- 1)^] 

Cone 

eo 

X 

VI 

X 

VI 

CM 

X 

[r? - (d - X) 2 ] 1/2 

irf-X 2 + 2dX + R 2 - d 2 ] 

Sphere 

X 3 < X < L 

R 1 

ttR 2 

Cylinder 


a Regions are defined by tangent point. 
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Input Instructions 


The input data for the computer code is entered in the following manner: 


Card 1 




Card column 

Variable 

Format 

Remarks 

1-72 

TITLE (I) 

(12A6) 

Any hollerith identification information 

Cards 2 , 5, 8, 

ofp 

V vv. 



1-4 

N 

(814) 

Number of intervals in Simpson's rule integration 
for determining X 7 . Must be odd number. 

Use N = 11 

5-8 

NN 


Number of time steps in later period. Use 
NN > 50 

9-12 

NNN 


Number of time steps in later period. Use 
NNN ^ 50 

Cards 3, 6, 9, 

etc. 



1-8 

R0 

(10E8.4) 

R q (see table I), ft 

9-16 

R1 


Rj (see table I), ft 

17-24 

D 


d (see table I), ft 

25-32 

XI 


Xj (see table I), ft 

33-40 

X2 


X 2 (see table I), ft 

41-48 

X3 


X 3 (see table I), ft 

49-56 

WP 


W p , flow rate, lb mass /sec 

57-64 

AL 


L, initial liquid level, ft 

65-72 

CONST 


Constant sometimes used to vary heating rates. 
Use 1.0 

73-80 

QW1 


2 

Average wall heat flux, Btu/(sec)(ft ) 
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Card 4, 7, 10, etc. 


Card column 

Variable 

Format 

Remarks 

1-8 

RHO 

(10E8. 4) 

p, density of propellant, tt) o /ft 3 

ill AO S 

9-16 

BETA 


Coefficient of thermal expansion of the propel- 
lant, V 1 

17-24 

AMU 


Viscosity of the propellant, lb m _ OC3 /(ft)(sec) 

IXiclSS 

25-32 

CP 


Cp, specific heat at constant pressure of the 
propellant, Btu/(°R)(lb mass ) 

33-40 

AK 


Thermal conductivity of propellant, 
Btu/ (sec) (°R) (ft) 

41-48 

TREF 


9 „, saturation temperature rise, °R 

fa 

49-56 

GLIL 


2 

g, acceleration due to gravity, ft/sec 

57-64 

Y 


Axial point in tank at which temperature history 
is desired, ft 

65-72 

XAMB 


Axial position of starting point of the boundary 
layer, ft 

73-80 

POWER 


Constant sometimes used for varying heating 
rates. Use 1.0 


As mentioned previously the integrated heating curves Q^Xg) and Q n (X g ) are also 
necessary input. Ordinarily, these curves are described in the subroutines WALL and 
NUC with tables of values of or Q n as a function of X g . The heating values are 
then determined from these tabular arrays by interpolation. However, at the option of 
the user, an equation or any other means may be used to describe the curves Q w and 
Q as a function of X . Communication with the main computer code is mainted through 
the two variables XXX and YYY which appear in common. The variable XXX corre- 
sponds to an axial coordinate, such as X c , and YYY corresponds to a heating rate, 

s 

such as Qw in WALL or Q n in NUC. The sample problem, presented in appendix B, 
gives one example of how these subroutines are used. 
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OUTPUT 


Variable 

Format 

Remarks 

AL 

(9E13.5) 

Same as input 

RO 


Same as input 

R1 


Same as input 

D 


Same as input 

XI 


Same as input 

X2 


Same as input 

X3 


Same as input 

WP 


Same as input 

GLIL 


Same as input 

RHO 

(9E13.5) 

Same as input 

BETA 


Same as input 

AMU 


Same as input 

CP 


Same as input 

AK 


Same as input 

TREF 


Same as input 

CONST 


Same as input 

QW1 


Same as input 

XAMB 


Same as input 

VI 

(8E13. 5) 

Volume in region of tank between 

0 < x< x v ft 3 

V2 


Volume in region of tank between 

Xj < x< x 2 , ft 3 
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Line 


4 


5,6,7, . 


Variable 

Format 

Remarks 

V3 


Volume in region of tank between 
X 2 < X s= X 3 , ft 3 

V4 


Volume in region of tank between 
X 3 < X< AL, ft 3 

T1 


Time required to outflow V^, sec 

T2 


Time required to outflow V 3 , sec 

T3 


Time required to outflow V 2 , sec 

T4 


Time required to outflow Vp sec 

PR 

(8E13.5) 

Prandtl number, dimensionless 

GT 


Grashof number divided by 

a 0-D ” 1 

V R 

TW 


V 0r 

H 


Heat transfer coefficient, Btu/ (ft 3 ) (sec) (°R) 

ETA 


n, dimensionless 

XSL 


x s,r ft 

POWER 


Same as input 

Y 


Same as input 

TIME 

(10E13.5) t, sec 

TEMP 


0(X,X 3 ), °R 

XS 


X s , ft 

F(XS) 


f(X„) or F(X„) for the initial or later period, 
s s 

respectively, dimensionless 

ACAP 


b(X s )/a(X g ) or y(X s )/a(X s ) for the initial 
or later period respectively, ft -1 
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Line 


Variable Format 


Remarks 


BCAP 

A < x s )Q w (X s ) _ 
c p w p e s a(x s ) 


A(X 1 > [<w + <w] 

CpVs “< X s> 


for the initial or later period respectively, 

. -1 
it ' 

QWAL 

Q w (X g ), Btu/sec 

QNVC 

Q n (X g ), Btu/sec 

VOLUME 

3 

volume of fluid in the tank, ft 

DEL 

6 (X g ), ft 


SAMPLE CALCULATION 

The sample problem described is for a 33 -foot-diameter nuclear -rocket propellant 
tank of liquid hydrogen with a geometry as shown in table I. The pressure head corres- 
ponds to the pressure difference of 22. 3 to 15. 0 psia. This is equivalent to a temperature 
rise of 2. 53° R or TREF = 2. 53° R. The flow rate is 3000 pounds mass per second. 

The heating rates, QNVC and QWAL correspond to nuclear radiation heat deposition in 
the propellant, and nuclear radiation heat deposition and ambient heating of the tank 
walls, respectively. 

The remaining input is presented in table II. The output for this calculation follows 
the computer listing in appendix B. A flow chart for the computer code is presented in 
figure 2. Note that the two subroutines, WALL and NUC, in the program listing are also 
input information. A plot of the resulting exit temperature history is presented in fig- 
ure 3. Typical running times average around 0. 04 minute per case on an IBM 709411/7044 
direct-couple system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 20, 1967, 

121-30-01-07-22. 
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TABLE H. - INPUT FOR SAMPLE PROBLEM 



NASA-C-836 (REV. 9-14*59) 













Calculate a(X s ), 
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Figure Z - Flow chart of computer code. 



























Temperature rise, 
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APPENDIX A 


SYMBOLS 


lgineering 

FORTRAN 

Definition 

symbol 

symbol 


a 

A 

variable defined by eq. (18) 

A 

AR, ARE, AREA 

2 

cross-sectional area of propellant tank, ft 

A' 

DAR, DARE, DAREA 

first derivative of A with respect to X g 

A” 

DDAR, DDARE 

second derivative of A with respect to X g 

A”' 


third derivative of A with respect to X g 

b 

B 

variable defined by eq. (17) 

C P 

CP 

specific heat at constant pressure, Btu/(°R)(lb magg 

d 

D 

see table I, ft 

F 

FXS 

variable defined by eq. (15) 

f 

FXS 

variable defined by eq. (2) 

i 


index used in summation in eq. (21) 

L 

AL 

initial liquid level, ft 

M 


upper limit defined in eq. (21) 

n 

AN, ETA 

parameter defined by eq. (6) 

Q 

QNVC 

integrated nuclear heating, Btu/sec 

Q w 

QWAL 

integrated wall heating, Btu/sec 

% 


o 

nuclear heating rate per unit volume, Btu/(sec)(ft ) 

% 


wall heat flux, Btu/(sec)(ft^) 

R 

TEST 

remainder defined by eq. (11) 

R 0 

RO 

see table I, ft 

R 1 

R1 

see table I, ft 

T 


temperature, °R 

T b 


bulk temperature, °R 

T i 


initial temperature of fluid, °R 
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saturation temperature corresponding to ullage pres- 
sure, °R 

TIME time, sec 

time when bottom of stratified layer first reaches X = 0, sec 
WP flow rate, lb mass /sec 

XI position of bottom of stratified layer, ft 

XI see table I, ft 

X2 see table I, ft 

X3 see table I, ft 

XS liquid level, ft 

XSL liquid level at t = tp, ft 
Y axial position in propellant, ft 

equivalent to left side of eq. (9) 
equivalent to right side of eq. (9) 

A variable defined by eq. (19) 

DEL stratified layer thickness, ft 

DELPR first derivative of 5 with respect to X c 

DELXS change in liquid level in one time inter Aral 
B variable defined by eq. (20) 

AMU viscosity, lb mass /( ft K sec ) 

defined by eq. (3) 

RHO density, lb mass /ft 3 

o 

surface area of tank, ft 
TEMP temperature rise (T - Tj), °R 

temperature rise of bulk fluid (T - Tj) , °R 
TREF saturation temperature rise (T g - T^), °R 

TW temperature rise across boundary layer, °R 
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APPENDIX B 


PROGRAM LISTING AND SAMPLE FOR TEMPERATURE PROFILES 


Computer Listing 


OIMENSICN XSFRUOU IOC) .TITLE 1 12) ,AR( 101) 

COMMON XXX, YYY 
lor FORMAT (814) 

102 fcrmat< icee.4 1 

110 FORMAT ( 1 HU 5X, 2H AL , 1 2 X * 2HR 0 , 1 1 X , 2H R l , 1 CX ♦ l HC * 1 2 X , 2 HX 1 , 1 1 X , 2 EX 2 * 1 IX 
1 , 2HX3 , 1 1 X , 2E WP , 10X,4HGLIL/(9E13.5) ) 

11 2 FTP mat ( 1H0,4X, 3HRHC, 1 1 X , 4H P FT A , 9X , 3HAMU. 10X.2HCF ,1 1X.2HAK, 11X,4ETR 
1 f F • 7X, 5FCCN5T * 1CX , 3H0W1 , 9X ,4HX AMP/ (9EI 3. 5 ) ) 

11 A FORMAT ( 1 HO * 5X f 2HV l , 1 2X , 2HV2 , 11 X , ?H V 3 * 1 C X ,2H V4 * 1 1 X , 2H T l , 1 IX , 2HT 2 ♦ 1 1 
1 X * 2HT3 , 1 IX , 2HT4/( 8E12.5I) 

Ilf FORMAT < 1FC,5X,2HPR,12X,2HGT ,1 1 > ,2HTV , 1 C X , 1HH f 1 2X ,3 LET A f 10X , 3HXSL *8 
lXt5HPDWERtl 1 X, 1HY/(8E13.5) ) 

11 8 FORMAT ( 1HO,4X,4HTIME,1CX,4HTEMF.10X,2HXS ,8X,5HF(XS),1CX,4HACAP.9X, 
1 4 FBC AP , 9X , 4HQWAL, 9 X , 4HC NVC , 8 X, 6H VC L LME ,1CX,3HDEL ) 

120 FORMAT ( 1CF12 .5 ) 

122 F0RMAT(12A6) 

RFAC <5, 122 ) (TITLECI) ,1 = 1.12) 

WRITF (6,122) (TITLEU >, 1 = 1,12) 

2 PF AC ( 5, ICC ) N , NN , KNN 

RfAO (5,102) R0,R1 ,D,X1 ,X2,X3, V*F , A L ,CC N S T ♦ CW 1 
R F AO (5,102) RHC, BETA, AMU, CF.AK, TREE, GLIL,Y,XAK8, POWER 
C 

C CALCULATE CONSTANTS FOR THIS RUN 

C 

KNN = C 
PI = 3.14159 
A 1 = -PI 
P T 2 = 2 .*P I 
PI = P I 2* R C 
A 2 = PI 

SC? = SORT (2. ) 

SC2 1 = $07-1. 

0 2 = R i*sc; 1 

0? = P I#RC*R0*$G21*SQ?1 
A3 = A 1 
pi = P 1 2*0 
P TR = PT*R 1*R1 
C 3 = P 1 R - P I *0*0 
C 4 s P I R 

PR = ( C P* AML ) / AK 

GT = RH0*RHC*GLIL*8ETA*AL*AL*At/(AML*AMU) 

T p M 2 = 0.13*AK/AL 
T F M 3 = (GT * P R ) * * 0 • 33333 
TW = (Ckl/(TEM2*TFM3))**C.75 
H = CVH/TVi 

AN = 4 • 0 C * ( 1RFF/TWI-1.C 
F 1 A - AN 
X 12 * X 1*X 1 
X 1 3 = X 1 2 * X 1 
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ATFC 


EFN SOURCE STATEMENT 


IFN(S) 


X22 * X2*X 2 
X 23 * X22*X2 
X 22 * X 3*X3 
X23 * >3 2*> 3 

VI = Al*X13/3.4Bl*X12/2. 

V2 = A2/3.*(X23-X13)4B2/2.*(X2 2-Xl2l4C2*<X2-Xn 
V3 = A3/3.*(X33-X23)4B3/2.*<X32-X22)4C2*<X3-X2 1 
V 4 - CA*(AL-X3) 

RWP * RHG/WP 
T 1 = VA*RWP 
T 2 = V 3 * R W P 
T 3 = V 2*R WP 
TA = V1*RWP 
A]N * A N 4 1 , 

A 2N = AN 4 2 • 

A 2N - AN43. 

AM * 1./A1N 
A fs 2 * AN1/A2N 
AN 3 * AN2/A 3N 

BFTAC = 1 • / (CP*WP*TREF*AN1 1 
XXX * XAMB 
CALL WALL 
OWX » YYY 
C 

C CALCULATE XSL 

C 

NXSL = 0 
H = N-l 
MN a M-l 

AINCR = 2 • *f AL-X3 ) 

X S PR ( 1 ) = X 3 

1C DFLXS * ( AL-XSPR ( 1 1) /FLOAT ( Ml 
AINCR = AINCR*. 5 
Of 11 TM,M 

11 X SPR C 1 4 1 J = X S P R { I 1 4 D E L X $ 

DC 12 1=1, N 

XXX = XSPRU) 

IFCXXX.GT.XA^B) GO TO 212 
OW( I) = O.C 
GC TO 12 
212 CALL WALL 

OWCI) = POWER* (YYY-OWX) 

1 2 CONTINUE 
IF(NXSL.GT.C) GO TO FC 

^ wi = oxu newm 

W 2 » 3.Q 
no 13 1=2, P . 2 

13 W 2 a W2 44 . C*QW ( I ) 

W 3 = O.C 

DC 1A I=3,RN,2 
1 A W 3 = W ?4 2 .C*OW ( I ) 

W = W14W24W3 

OFL = RFTAC4(DELXS/3.C )*W 
IFCNXSL.G7.C ) GO TO 
TF ST = OEl-XSPR(l) 

76 IF (AB5( IFST)“.CC1I 18,18,15 
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ATHC - EFN SOURCE STATEMENT - IFN(S) - 

1* KNN = KNK+1 

I F ( KNN *GE • 5 C ) GO TO 25 
I F ( TEST 116,18,11 

16 IF<KNN.FC. 1 1 GO TO El 

IF( INXSL.FC.il . AND • ( KNN .EQ. 2 ) 1 GO TC 82 
IF( (NXSL.EC.2) .AND. (KNN.FG.3) ) GO TC 83 
X5PRU) = X SPR ( 1 1-C.5+AINCR 
GC TO 10 

17 XSPRC1I * XSPR(1)+C.5*AINCR 
GC TO 10 

25 X 51 * XSPR(l) 

GC TO 26 

PI NXSL = l 

XFPR<1) = X2 
AINCR = 2 • * ( X3-X 2 ) 

GC TO 10 

82 N > SL = 2 
XSPR(l) = >1 
AINCR = 2 • * ( X 2- XI ) 

GC TO 10 

83 NXSL ■= 3 

X 5PR ( 1 ) = C.O 
AINCR = 2 • ♦ X 1 
GT TO 10 
RC nc 84 1*1, N 
7 = XSPR t I 1 
IF(7.GE.X3)GO TO 85 
IF(7.CE.X21GO TO 86 
I F( 7 .GF • X 1 ) GO TO 87 
AR< H = A 1* 7^7 !♦ 7 

GC TO 88 

87 A P ( T 1 * A2*Z*Z + B2*‘Z + C2 
GC TO 88 

86 A P ( I 1 = A 3* Z^ Z+B 3^ Z+C 3 
GC TO 88 

85 A R ( I ) = C4 

8 P CONTINUE 

owm * cwm*ARU i 

84 CONTINUE 
GC TO 89 

899 IFtXSPRm.GE.X3J GO TC 95 
I F ( XSPR ( 1 1 . GE • X2 1 GO TO 96 
IF(XSPRU).GE.Xl) GO TO 97 
OAR = P I 2* ( -XS PR ( 1 > ♦RO ) 

OTAR = -PI2 
GO TO 99 

97 CAR = PI2* ( XSPR ( 1 ) +RC*SC21> 

OCAR = P 1 2 
GC TO 99 

96 DAP = PI2*(-XSPR(1 1*0) 

DCAR = -PT2 
GO TO 59 

95 car = c.o 
c r ap = c.o 

99 7 7 = (X SP P ( 1 J *A R ( 1 ) * AN 1 -X SPR (1) ♦XSPR <11 ♦AK2*CAR*X SP P < 1 1 ♦XSPR <11 ♦XSP 
1 R ( I ) ♦DC AP^ AN3 ) /AN1 
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nnn n r> 


ATHC 


FFN SCURCE STATEMENT 


IFN(S) 


TEST = DEl-ZZ 
GC TO 76 
IP X5L = XSPR(l) 

r 

X SL CALCULATION FINISHED 

7 6 7 FT A 1 = AN1*XSL 
Z ETA2=AN2 * XSL+XSL 
7FTA3 = AN3*XSL*X$L*X$L 
ALPHAN = F I 2*7 ETA3 
BETAN = FI2*ZETA2 
GAMMAN = FI2*R0*ZETA1 
E TAN = RC*B ETAN 
DEL = C.C 

DFLXS = (Al-XSL) /FLOAT (NN) 

WRITF (6,110) AL,RC*R1,D,X1»X2»X3,WP,GLIL 
WRITE (6,112) RHO, BETA, AMU, CP, AK,TREF, CONST, QWl,)(ANe 
WRITE (6,11M Vl,V2,V3,V4,Tl,T2,T3,T4 
WRITE (6,116) PR, GT,TW,H,FTA,X$L, POWER, Y 
WRITE (6,1 IB) 

XS = AL 
GO TO IS 

201 0 EL PR = -BETAO+OWAL 
FXS * C.O 
OCWAL = C.O 
N7 2 0 

A = V 1 +V 2 +V3+V4 
B AN 1 = BET AC*AN l 
B = -QWAL*C4*BAN1 
ACAP = B/A 

BCAP =C4*CNVC*BANl/A 
20 XS = XS-DELXS 
1 9 XXX = XS 

CALCULATE HFATING RATES 
CALL WALL 

IF (XS-XAMB) 71,71,72 

71 ONX = YYY 
OWX = C.C 
GO TO 13 

72 QM = C.C 
OWAL - PC WER* ( YYY-QWX I 
GC TO 77 

73 C W AL = O.C 

77 CALL NUC 

QNVC = POVER*(YYY*CONST+CNX+CWX) 

XS? = XS+XS 
X S3 = XS2+XS 

I F ( XS .GT • AL-OFLXS/2. ) GC TO 201 
IFCXS.GT.X2 ) GO TO 9C 
IF(XS.GT.X2) GO TO 91 
TFIXS.GT.Xl ) GO TO 92 
ARE = A1*XS?+B1*XS 
CARE = PI2*(-XS+R0) 

C T AR E = -P 12 
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n rt n n n r> r> 


ATHC 


PFN SOURCE STATEMENT 


IFNI S) 


OMEGA * A l*XS3/3.+Bl*X$2/2. 

T IMF = T14T2 + T3+RWPMV1-CMFGA) 

GO TO 94 

9 2 ARE * A2*XS2+B2*XS+C2 
CARE = PI2*(XS+RO*SC21 ) 

C T ARE * P I 2 

CMFGA = V 1 * A2* ( X5 3- X 1 3 ) /3 • ♦ B2^ (XS2-X12 ) / 2 • +C2* (XS-X1 ) 

TIME = T14T2+RWP*(V1+V2-0MEGA) 

GC TO 94 

91 ARE = A3*XS2+R3*XS+C3 
H AR F = P I 2* (— X5+D ) 

CCARF = -P 1 2 

OMEGA a V1 + V2 + A3M XS3-X23) /3. + B3M XS 2- X2 2 ) /2 . +C3 ♦ ( XS-X2) 

TIME = THPWP*( V1 + V2 + N/3-CMEGA) 

GC TO 94 
9 C ARE = C4 
CARE = C.C 
CCARF * 0,C 

OMEGA = V 1 +V 2+ V 3+C 4* ( X 5-X3 ) 

TIME = C4*RH0*( AL-XS » / %r P 
94 AREA = ARE 
CARFA * CARE 
IF (XSL-XS) 2? • 2 2 * 46 

CALCULATE INITIAL PERIOO 

22 CONTINUE 

0 EL = DEL-OELPR*DELXS 
DEL? a OEl*DEL 

DEL PR a - ( ARF*BAN1*QWAL+0EL*CAFE 3 *AN1-DEL2*CCARE*AN2 ) /( ARE*ANl-2 
1EL«rARE^AN243.^0EL2»CnARF*AN3) 

I F(CEL •€ E • X S IDELaXS 
XI = XS-OFL 
PP = B 

B = -ARFA*CWAL*BAN1 

OCWAt = COWAL-(B>BP)*OELXS/2. 

A = CMFGA-CQWAL 

3 P 0 = (AREA*CNVC )/(CP*WP*TREF ) 

FXS a FXS+(BCAP-ACAP*FXS)*DELXS 
ACAP = 6/A 
BCAP = O/A 
IF (Y-XI ) 40 *40*42 
40 TEMP = TREF*FX S 
GC TO 44 

42 TFM1 = ( ( Y-Xl )/OEL)**AN 

TEMP = TREF*(FXS*(KC-TEM1 ) 4 ( TEM1 ) ) 

44 I F( TFMP.GF.TRFF >TEMP = TRFF 

WRITE (6*120) TTMFfTEMP*XStFXStACAPtBCAP*OWAL*CNVCtOMEGA*OEL 
GC TO 20 

INITIAL PERIOO FINISHED 
CALCULATE FINAL PERIOD 

46 PHI = ( 1 .-XS/XSL l**AN 
IF(XS*LE.Xl) GO TO 48 
I F (XS.LF .X2 ) GO TO *0 
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ATHC - EFN SOURCE STATEMENT - IFN(S) - 

IF<XS.LF.X3) GO TO 52 
XS^X = (XS-X3 1 /XSL 
E 33 = C ] .-XS3X)**A3N 
F 23 = (1 .-XS3X)**A2N 

52 X S2X = (XS-X21/XSL 
E32 = < 1 .-XS2X)**A3N 
F 2 2 = ( 1 .-XS2X )**A2N 

5C XS1X = (XS-Xll/XSl 
F31 = ( 1.-XS1X)**A3N 
F 21 = (1 XS IX )** A 2N 

AS XSX = XS /XSL 

E ’C = (1 .-XSX)**A3N 
F2C = ( 1 .-X$X)**A2N 
F 10 = < 1 .-XSX)**A1N 
IF(XS.GT.X3) GO TO 54 
TF(XS.GT.X2 ) GO TC 56 
TFtXS.GT.Xl ) GO TO 58 
IFtXS.GT.C.) GO TO 6C 
GO TO 2 

6C A = 0MEGA-7ETA1*AREA47FTA2*CAREA— ALPHAN*(-1.0+F30>-FTAN*E2C 
B = 7ETA1*CAREA4-BETAN*(1.C-E20I-GAMMAN*E10-PHI 4AREA 
GO TO 62 

58 A = 0MEGA-7ETA1*AREA+7ETA2*CAR E A— A L PHAN* (1.0-2.0*E?1+E30)-ETAN*E20 
B = 7ETA1*CARFA+BETAN*C-1.C+2.C*E21-E?0) -GAM MAN* ElO-PHI+AREA 
GO TO 62 

56 A = OMEGA-7ETA1*AREA47ETA2*CAREA-ALPHAN*(-1.0+2.0*E32-2.0*E31 4 
1 F 30 ) - ETAN*E20 

B = 7FTA1*CAREA 4 EETAN*< 1 . C-2 .C*F2242 .0*E21-E2C > -GAMM AN*E 10- 
1PFI+AREA 
GC TO 62 

54 A = OMFGA-2ETA1*AREA-ALPHAN*(-E3342.0*E32-2.0*E314-E30)-ETAN*E20 
B = BFTAN*(F23-2.C*E2242.0*E21-E20)-GAMMAN*E10-PHI *ARFA 

62 0 = ( AREA*(CNVC4QWAL >1 / (CP*V>P*TREF 1 
mN7.GT.01G0 TO 2 CC 
ACAP = B/A 
BfAP = 0/A 
N 7 = 1 

2QC F XS = FXS4(BCAP4-ACAP*( l.C-FXS) )*OFLXS 
ACAP = p/A 
BCAP = 0/A 
TF (V-XS1 64.64,66 

64 TEMP = TREF*(FXS*(1.C-(1.0-(XS-Y)/XSL)**AN)4(1.C-(X S-Y l/XSL ) **AN» 
GO TO 68 

66 TEMP = TREF 

68 1FCTFMP.CF.TREF1TFMP = TRFF 

HP I TF (6,1201 TIMF, TEMP, XS.FXS, ACAP, BCAP .OWAL.CAVC, OMEGA 
CELXS = X S L/FLOAT ( NNM 
GC TO 20 
ENC 
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S IPROUT I NE WALL 
CCMMCN XXX, YYY 

YYY * 2.565 89*(XXX/33.C- 27. C/64. C) 
1F1XXX.LE.13.9219) YYY = C.C 
PFTURN 
END 


SUBROUTINE NUC 
DIMENSION XI13I.VC13* 

DATA X,Y/C.C*. 66 667, 1.47,2.3, 5. 6, 8. 9,12.1, 15. 6,18.9,22. 3, 25.6*68. 9 
1 , 100.0,0.0,4 2. 68 3, 13.415,102.440,215.123,312.441,384.149,4 28.5 39,4 
241. 320, 457.564, 462. 6F6 ,465.515 ,469.515/ 

CCWMON XXX, YYY 

cn ic 1=1,13 

TF(XXX.LE.X( II) GO TC 20 
1C CCNTINUE 
2C M = 1-1 

I F ( M )4C,4C ,30 

3C YYY = Y<M + (Y(M+1 ) - Y C « )) * ( XXX- X ( H 1 } / { X( M+l )-X m ) 

GO TO 50 
40 YYY = C.C 
50 RETURN 
END 
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03 
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03 
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C3 

0, 
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RO 

C.45752F 01 
BETA 

C • 10000F-C 1 
V 2 

C.16167E C4 


OT 

C .4C10 IE 17 
TEMP 

0*159086-02 
0.31857E-02 
C.47848E-C2 
C.63881F-C2 
C.75956E-02 
C.96074E-C2 
0,1 1223E-C1 
C.12844E-01 
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0.243C8F-CI 
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C2 
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C 5 

V3 


0.83153F 

CA 
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cc 
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(2 
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C 2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 

C.65178E 

C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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C2 
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(2 
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C2 
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C2 
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C2 
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C2 

C.63C P9E 

C2 

C • 62 929F 

C 2 


D 
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C 2 

CP 
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Cl 

VA 
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C5 

H 
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-C 2 


Sample 

XI 
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AK 

0 • 35000 E— 05 
T 1 
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ETA 

C.33266F C2 


F ( XS I 
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0. 960786-02 
0 • 10262E-C l 
C.10918E-01 
0. 115766-01 
0* 12236E-G1 
C.12897E-01 
C. 13561E-01 
0*142266-01 
C.1A893E-01 
C • 1556 IE— 0 1 
C.16232E-01 
C.169C5E-01 
C.17579E-C1 
0.18255E-01 
0. 18934E-01 
C* 19614E-01 
C.20256E-01 
0# 2098 OE -01 
C.21666E-01 
0*22 35 36-01 
C • 2 3043E-0 1 
C. 237356-01 
C.24429E-01 
0. 251256-01 
C. 250226-01 
0.26522E-01 
0.27224E-01 
0.279286-01 
0.28634E-01 
C. 293426-01 
C.30052E-01 
0* 30764 F-0 1 
0 • 3 147 8E — 0 1 
C.32155F-01 
0.32913E-01 
0.33634E-01 
0.34357E-01 
0* 350826-01 
C.358C9E-01 
C * 36 5 38F -01 


ACAP 

-0*372296-04 
0 *37 322E-04 
C * 373 1 A F-CA 
-0*373C7E-04 
-C *37300 E-OA 
-0 • 37292E-04 
-0 *3728 5E-G4 
-0.37278E-04 
-0.37270E-04 
-0.37263E-04 
-0. 372556-04 
-0*37248 E-OA 
-0 *37240 E-OA 
-0 *37232 E-OA 
-0 *372 25E-04 
-0 • 372 1 7 E-OA 
-C . 37 2C9E-04 
-0 . 3720 1E-0A 
-0 • 37 1 93E-04 
-0.37185E-CA 
-0.371 77E-0A 
-0 • 37169 E-OA 
-0*37161 E-OA 
-0 • 371 53E-0A 
-0 *37 1 45 E-OA 
-0,371376-04 
-0 ,371296-04 
-C *37120 E-OA 
-0*371 1 2E-0 A 
-0*371 CAE-OA 
-0.37CS5E-04 
-C.37C87E-04 
-0.3 70 78E-04 
-0.37C70E-CA 
-C.37C61E-0A 
-0.37052E-CA 
-C.370A3E-CA 
-C.37035E-0A 
- C • 370 26 E— 04 
-0 *3 7017 E-OA 
-C.370C8E-0A 
-0 • 36999E-0 A 
-0 • 36990 E-OA 
-0.36981 E-OA 
-0.36572E-C4 
-0 # 3696 3 E-C A 
-0 • 369 5 3E-0A 
-0 • 369 A A E-C A 
-C # 36935E-04 
-0*3 6925 E-OA 
-0,369166-04 
-0 • 36906 F-OA 
-0 • 3 6857E-C4 
-0 • 36887 E-CA 


PRECEDING PAGE BLANK NOT RLJVu— < 


Output 


X2 

X 3 


V»P 


GLIL 




C .96C30E 01 

0.21269E 

02 

C.3CCC0E 

C3 

C.32200E 

02 



TREE 

CONST 


owi 


XAMB 




0 • 25 300E 01 

0.1C0C0E 

Cl 

C • 7 5CC0E-C3 

C.140CCE 

02 



T2 

T3 


T4 






0.11946E 03 

0.23226E 

02 

C.43132E 

CO 





XSL 

POWER 


Y 






0.73434E 01 

0.1C000E 

Cl 

C. 






0CAP 

OWAL 


CNVC 


VOLUME 

OEL 

0.3924 IE-02 

0.44665E 

01 

C.46952E 

C3 

C.52876E 

C5 

0. 12739E-01 

0.3934 3E-02 

0.44540E 

01 

0.46952E 

C3 

C.52739E 

05 

0. 25443E-01 

0.39446E-02 

0.44415E 

Cl 

0.46952E 

C3 

C.526C1E 

C5 

0.3811 IE- 

01 

0.39550E-02 

0.44290E 

Cl 

C.46S52E 

C3 

0.52464E 

05 

0. 50744E-01 

0.39654E-02 

0.44165E 

01 

C.46952E 

C3 

C .523266 

05 

0.63341E- 

•Cl 

0 • 3975 8E-Q2 

0.44040E 

01 

0.46S52E 

03 

0.52189E 

05 

0.75902E- 

01 

0.39 863E-02 

0.43915E 

01 

C.46952E 

C3 

C.52052E 

05 

0.88428E-01 

0.39969E-02 

0.43790E 

01 

C.46S52E 

03 

0.51914E 

05 

0. 10092 E 

00 

0.40076E-02 

0.43665E 

01 

0.46952E 

C3 

C.51777E 

05 

0.11337E 

00 

0.40182E-02 

0.43540E 

01 

0.46952E 

C3 

C.51639E 

05 

0.12579E 

00 

0.4C290E-0 2 

0.43416E 

01 

C.46952E 

C3 

C.515C2E 

05 

0.13818E 

00 

0.40398E-02 

0.43291E 

01 

C.46952E 

C3 

C.51365E 

05 

0.15053E 

00 

0. 4050 7E— 02 

0.43166E 

01 

C.46952E 

C3 

C.51227E 

05 

0.16284E 

00 

0 .406 16E-02 

0.43041E 

01 

C.46952E 

C3 

C.5109CE 

05 

0.17512E 

00 

0.4C725E-02 

0.42916E 

01 

C.46952E 

C3 

C.5C952E 

05 

0.18736E 

00 

0.40836E-02 

0.42791E 

01 

C.46952E 

03 

C.50815E 

05 

C.19956E 

CO 

0.40947E-02 

0 .4 26666 

01 

C.46952E 

C3 

C.50678E 

05 

0.21173E 

00 

0.41058E-02 

0 .4254 IE 

Cl 

0.46952E 

C3 

0.5C540E 

05 

0 .2238 7 E 

00 

0 .4 1 17 IE- 02 

0.424 16E 

01 

C.46952E 

C3 

C.504C3E 

05 

0.23597E 

00 

0.4128 3E-02 

0.42291E 

01 

0.46952E 

03 

0.5C265E 

05 

0.24803E 

00 

0.4 1397E— 02 

0.42166E 

01 

C.46952E 

C3 

0 .501286 

05 

C.26006E 

00 

0.41511 E— 02 

0.42041E 

Cl 

C.46952E 

C3 

C.4999CE 

05 

0.27206E 

00 

0 .4 1626E-02 

0.41917E 

Cl 

C.4E952E 

C3 

C.49853E 

05 

0.28401E 

00 

0 .4 174 IE-02 

0.41792E 

Cl 

C.46952E 

C3 

C.49716E 

05 

0.29593E 

00 

0.4 1857E-C2 

0.41667E 

01 

C.46952E 

C3 

C.49578E 

05 

0.30782E 

00 

0*4197 3E-02 

0.41542E 

Cl 

C.46952E 

C3 

C.49441E 

05 

0.31967E 

00 

0.4209 IE-02 

0.41417E 

01 

C.46952E 

03 

C.49303E 

05 

0.33149E 

00 

0.4220 9 E— 02 

0.41292E 

Cl 

C.46952E 

C3 

C.49166E 

05 

0.34327E 

CO 

7C C2 

£ _ A 1 1 A7P 

n i 

0.46S52E 

C3 

C.49C29E 

05 

0.25501E 

00 

0.42446E-02 

0.41042E 

01 

C.46952E 

C3 

C.48891E 

05 

0. 36 672 E 

00 

0.42566E-02 

0 .40917 E 

Cl 

C .46952E 

C3 

C.48754E 

05 

0.37839E 

00 

0.4 268 7E-02 

0.4C792E 

01 

C.46952E 

C3 

0.48616E 

05 

0.39003E 

00 

0.42808E-0 2 

0.4C667E 

Cl 

C.46952E 

C3 

C.48479E 

05 

0.40163E 

00 

0 .4 293 0E-02 

0.40542E 

01 

C.46952E 

C3 

C.48342E 

05 

0.41320E 

00 

0.4305 3 E— 02 

0.4C418E 

Cl 

C.46952E 

C3 

C.482C4E 

05 

0.42473E 

00 

0.43 176E-02 

0.4C293E 

Cl 

C.46952E 

C3 

C.48067E 

05 

0.43623E 

00 

0.43300E-02 

0.4C168E 

01 

C.46952E 

C3 

0.47929E 

05 

0.44769E 

00 

0 .4342 5E-02 

C.40043E 

01 

C.46952E 

03 

0.47792E 

05 

0.4591 1 E 

00 

0.43550E-02 

0.39918E 

01 

C.46952E 

03 

C.47655E 

05 

0.47050E 

CO 

0. 43677E-02 

0.39793E 

01 

C.46952E 

C3 

C.47517E 

05 

0. 48186E 

00 

0 .4 380 4E-02 

0.39668E 

01 

0.46952E 

C3 

C.47380E 

05 

0 • 49 3 1 7 E 

00 

0.4393 1 E-02 

0.39543E 

01 

C.46952E 

C 3 

C.47242E 

05 

0.50446E 

00 

C.44C60E-02 

0.39418E 

01 

C.46952E 

03 

C.471C5E 

05 

0.5157CE 

00 

0.44 189E-02 

0.39293E 

01 

C.46952E 

C3 

C.46968E 

05 

0.52691E 

00 

0.4431 9E-0 2 

0 • 391 68E 

Cl 

C.46952E 

C3 

0.46830E 

05 

0.53809E 

CO 

0.4445 0 E-02 

0.39044E 

Cl 

C.46S52E 

C3 

C.46693E 

05 

0.54923E 

00 

0.4458 1 E-02 

0.38919E 

01 

C.46S52E 

C3 

C.46555E 

05 

0.56034E 

00 

0 .4471 3E-02 

0 . 387S4E 

Cl 

C.4E952E 

C3 

C.46418E 

05 

0.57140E 

00 

0 • 44846E-02 

0.38669E 

01 

0.4*9526 

C3 

C. 46281 E 

05 

0.58244E 

cc 

0 .4498 0E-02 

0.38544E 

Cl 

0.46952E 

03 

0.46143E 

05 

0. 59 344 E 

00 

0.4511 5E-02 

0.38419E 

Cl 

C.46952E 

C3 

C.46C06E 

05 

0.60440E 

CO 

0.45250E-02 

0.38294E 

01 

0.46S52E 

C3 

0.45868E 

05 

0.61533E 

00 

0.4538 7E-02 

0.38169E 

01 

C.46952E 

C3 

C.45731E 

05 

0.62622E 

00 

0.45524E-02 

0.38044E 

Cl 

C.46952E 

03 

C.45594E 

05 

0.63707E 

00 


29 



AL 


0.716C4E 

C? 

RHO 


0.4310CE 

Cl 

VI 


C.3CC22E 

C2 

PR 


0.5391 4E 

Oi 

TIME 


0.1C857E 

C3 

C.l 1055E 

03 

0.11252F 

C3 

0.1 145CF 

C3 

0.1 16476 

C 3 

0 • 1 184 4 E 

C 3 

C.l 204 2F 

C3 

0.12239F 

03 

0.1 243 76 

C3 

0.126346 

03 

0. 128316 

03 

C. 120296 

C3 

0.132266 

03 

0.13424F 

03 

0 . 1 362 1 6 

03 

0.13818F 

03 

C. 140166 

C 3 

0.142136 

03 

0.1441 IF 

03 

0.146C6F 

C3 

0.148056 

C 3 

C. 150036 

03 

0.1520C6 

03 

0.153986 

03 

0.155956 

C 3 

0.15792F 

03 

C.159SCE 

03 

0. 161876 

C 3 

0.16385F 

03 

0.165826 

03 

0.16775E 

C3 

0.169776 

C 3 

C. 17174E 

C 3 

0.17372E 

03 

0.1756SF 

03 

0.17766E 

03 

C. 17964F 

03 

0.181616 

03 

0. 1P359E 

03 

0. 1 855 6 F 

C3 

0.187546 

C3 

0.1 895 1 E 

C3 

0.191486 

C3 

0.193466 

03 

0.195436 

C3 

0. 197416 

C3 

C. 199386 

C 3 

0.201356 

C 3 

0. 20333E 

C3 

0. 2053 OF 

Cl 

0.207286 

03 

0. 209256 

03 

0.2 11226 

03 


RO 

C.49792E 01 
BETA 

C.10090F-C1 


V 2 


C. 161676 

C4 

C» T 


C.4C131E 

1 7 

T c MP 


0.94292E- 

■Cl 

C. 961496- 

-Cl 

C .980116- 

-01 

C.99379E- 

-Cl 

0.1C175E 

CO 

0. 103636 

cc 

0. 10552E 

00 

0. 1C 7 416 

cc 

0. 1 093 OF 

00 

0.111216 

00 

C. 11312E 

CO 

C.11503E 

GO 

0.11695E 

00 

C. 118886 

oc 

C. 120816 

00 

0. 122756 

00 

0. 124696 

oc 

C. 1266 56 

00 

0.1286C6 

00 

0.13057E 

CO 

C.13254E 

CO 

0. 134516 

00 

0.13650E 

00 

0.13848E 

CO 

0. 14048E 

00 

C. 14248E 

cc 

C.14449E 

cc 

C. 1 46 5 IE 

cc 

0. 14853E 

00 

C.15056E 

CO 

0. 15259E 

oc 

0.15464E 

oc 

C.l 566 9E 

00 

C.15874E 

cc 

0. 160806 

cc 

0. 16287E 

00 

C. 164956 

oc 

C.16704F 

00 

C.16913E 

00 

C.17122E 

CO 

0. 17333E 

00 

0. 17544E 

CO 

0.17756E 

00 

C.17969E 

00 

C. 181836 

00 

C. 16397E 

oc 

C. 18612E 

cc 

C.18828F 

CO 

0. 19044F 

oc 

0. 19261E 

oc 

C.15479E 

cc 

C. 196986 

oc 

C. 19918E 

CO 


R 1 


C • 1650 CE 

C2 

AMU 


C.7400GE- 

C5 

V3 


0.83153E 

C4 

TW 


C.29534E 

CC 

XS 


0.627686 

C2 

C.626C7E 

C2 

0.62447E 

C2 

C • 622866 

C2 

C.62126E 

C2 

C.61965E 

C2 

C.618C4E 

C2 

0.61644E 

C2 

C. 614836 

C2 

0.613226 

C2 

C. 611626 

C2 

C.610C1E 

C2 

C.6C04OE 

C2 

0.60680E 

C2 

0.605196 

C2 

0.603586 

C2 

C.6C19BE 

02 

C.60037E 

(2 

C. 598766 

C2 

C.59716E 

C2 

0.595556 

C2 

0.593946 

C2 

C.59234E 

C2 

C. 590736 

C2 

0.589126 

C2 

C. 587526 

C2 

0.58591E 

C2 

0 • 5843 1 E 

C2 

0.582706 

C2 

C.581C9E 

C2 

C.57949E 

C2 

C .577886 

C2 

C.57627E 

C2 

0 .574676 

C2 

C.573C6E 

(2 

0.571456 

C2 

0.56985E 

C? 

C.56824E 

C2 

0.566636 

C2 

0.565036 

C2 

C. 563426 

C2 

C. 561816 

C2 

C. 560216 

02 

C.5566CE 

C2 

C. 556996 

C2 

C.55539E 

C2 

C.55378E 

C2 

0.55218E 

C2 

C. 550576 

C? 

C. 548966 

C2 

0.54736E 

C2 

0.54575E 

C2 

0.54414E 

C 2 


0 

C. 212696 C2 
CP 

C.25500E Cl 
V4 

0 • 43C5 IE C5 
H 

C.25395E-C2 

F(XS> 

C.37270E-01 
C • 380 C4E-01 
C .387406-01 
C • 394 7 6E -Cl 
C.40218E-01 
C .4096 1E-0 i 
0.41 7 C66-C 1 
0.42454E-01 
C.43203E-01 
0 .4 3955E-01 
C.44710E-01 
C.45466E-01 
Q.46226E-01 
0.46987E-0 1 
0.47 7516 -01 
C. 485176-01 
C. 492866-01 
0.500576-01 
0.50831E-01 
0.5 16C7E -0 1 
0. 523866-01 
0.531676-01 
0. 539516-01 
C. 54737E-01 
C.55526E-01 
C. 563 17E-01 
0.571116-01 
0.579C8E-01 
C . 587 C76-0 1 
C.595C9E-01 
0.603146-01 
C. 61 1216-01 
0.61931E-01 
0.62744E-01 
C . 63559E-0 1 
0.643776-01 
0.65 198E-0 1 
C.66022E-01 
C.66848E-01 
0.676 78E-0 1 
C. 685106-01 
C • 69345E-0 1 
0.70183E-01 
C. 710246-01 
0.710686-01 
C.727 15E-C1 
C. 735646-01 
0.744176-01 
C. 752736-01 
0.761326-01 
C. 769946-01 
0.778596-01 
0 .787276-01 


XI 

0.14584E 01 
AK 

0 • 35000 E— 05 
T 1 

C. 618516 03 


ETA 

C.32266E C2 
ACAP 

-0.36877E-C4 
-0 .36867 E-C 4 
-0.3685 8 E-C4 
-0.36848E-G4 
-C.36838 E-C4 
-0.36828E-04 
-0.36818E-04 
-0.368C7E-04 
-0.36797E-04 
-0 • 36787E-04 
-0.36777E-C4 
-0.36766E-04 
-0.36756E-C4 
-0.36745E-C4 
-0.367356-04 
-0.367 24 E-04 
-0.36713 E-04 
-0 .367026-04 
-0.36691E-C4 
-0 .36680 E-04 
-0.36669E-04 
-0.366586-C4 
-0.366476-04 
-0 • 366 36E-C4 
-0.366256-04 
-0.36613E-C4 
-0.36602E-C4 
-0.365906-04 
-0 • 365 796-C4 
-0.36567.E-04 
-0 . 36555E-C4 
-0.365436-04 
-0 .36531 E-04 
-0.365 19E-04 
-0 • 36507 E-C4 
-0.36495E-04 
-0.36483E-C4 
-0.36470 6-C4 
-0.36458E-C4 
-0.36445E-C4 
-0 • 364 33 E-04 
-0 .3 64 20E-C4 
-0 .36407 E— C4 
-0.363946-04 
-0 • 36 38 1 E-04 
-0.36 368 E-04 
-0.36355E-C4 
-0 • 36 342F- C4 
-0.3 6 3 29 E-04 
-C *363 1 5 E-04 
-0.363C2E-04 
-0 • 36288 E-04 
-0.36274E-C4 


X2 


C.96C30E 

01 

TREF 
0.25 300E 

01 

T2 

0.11946E 

03 

XSL 

0.734346 

01 

BCAP 



C .456628-02 
C.45801E-02 
0.45940E-02 
C .4608 IE— 02 
C.46222E-02 
C.46364E-02 
C.46507E-02 
0. 46 651 E- 02 
0 .46796E-02 
0.4694 2E— 02 
0.47C88E-02 
0.47236E-02 
0.47384E-02 
C.47534E-02 
0.47 68 4E-02 
0.47836E-C2 
C.47S88E— 02 
0 .48 141 E-02 
C.48295E-02 
0.48451 E-02 
0.4860 7E-02 
0.46764E-02 
0.4892 2E-0 2 
C.49082E-02 
C • 4924 2E-02 
0.49404E-02 
0.49566E-02 
C . 49 730E-02 
0.49894E-02 
C • 5006 C E-02 
C . 5022 7E-02 
0 • 5039 5E-02 
C.5C564E-C2 
0.5C734E-02 
0 • 5090 5E-02 
0. 51C78E-02 
0.51251E-02 
0 . 51426E-02 
0 . 5 160 2E-02 
0.5 1779E-02 
0.51958E-C2 
C.52138E-02 
0 . 5 231 8E-02 
0 .5250 IE-02 
0.5 2684E-C 2 
0 . 5 286 9E-0 2 
0. 53055E-02 
C.53242E-02 
0 .5343 1 E-02 
0 .5362 1 E-02 
C.53813E-02 
0.54C05E-02 
0.54199E-C2 


X3 


0.21269E 

02 

CONST 


O.ICOCOE 

Cl 

T3 


0.22226E 

02 

POWER 


O.ICOCOE 

Cl 

QWAL 


0.37919E 

01 

0.37794E 

Cl 

0.37669E 

01 

0 .375456 

01 

C.37420E 

Cl 

0 • 372S5E 

01 

0.37170E 

01 

C.37045E 

01 

0.36920E 

01 

0.36795E 

01 

0.36670E 

Cl 

0.36545E 

01 

0.36420E 

01 

C.36295E 

01 

0. 36171 £ 

01 

C.36046E 

01 

0.35921E 

01 

0 • 35796 E 

01 

0.35671E 

01 

0.35546E 

01 

0.35421E 

Cl 

0.35256E 

Cl 

0.35171E 

Cl 

0.35046E 

Cl 

0.34921E 

Cl 

0.34756E 

Cl 

0.34672E 

Cl 

0.34547E 

Cl 

0.34422E 

01 

0.34297E 

01 

0.34172E 

Cl 

0.34047E 

01 

0.33522E 

Cl 

0.33757E 

Cl 

0 .3 3672E 

Cl 

0.33547E 

01 

0.33422E 

Cl 

0.33297E 

01 

0 . 33 1 73 E 

Cl 

0.33048E 

01 

0 .329236 

Cl 

0.32798E 

Cl 

0.32673E 

Cl 

0.32548E 

Cl 

0.32423E 

01 

0 - 32298E 

01 

0.32173E 

01 

0.32048E 

01 

0.31923E 

01 

0.31799E 

Cl 

0.3 1674E 

Cl 

0.31549E 

Cl 

0.31424E 

Cl 


WP 

C.3CCC0E C3 
OKI' 

C .7 5CC0E-C3 
74 

C.42132E CO 


Y 

C. 


CNVC 


0.46952E 

C3 

C.46952E 

03 

G.46952E 

C3 

0.46952E 

C3 

C.46952E 

03 

C.46952E 

03 

C .46952E 

C3 

C.46952E 

03 

C.46952E 

03 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

C3 

C.46S52E 

C3 

C.46952E 

03 

0.46952E 

03 

C.46S52E 

C3 

C.46S52E 

C3 

C.46952E 

03 

C.46952E 

C3 

C.46952E 

C3 

C.46S52E 

03 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

03 

C.46952E 

C3 

C.46952E 

03 

0.46952E 

C3 

C.46952E 

C3 

C.46952E 

03 

C.46952E 

C3 

0.46952E 

C3 

C.46952E 

C3 

C.46952E 

03 

C.46952E 

C3 

0.46952E 

03 

C.46952E 

C3 

C .46952E 

03 

C.46952E 

C3 

G.46952E 

C3 

C.46952E 

C3 

C.46952E 

03 

C.46952E 

C3 

C.46952E 

C3 

C.46S52E 

C3 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

C3 

C.46952E 

C3 

0.46952E 

03 

C .469 52 E 

C3 

C.46952E 

C3 


GLIL 

C.32200E 02 

XAMB 

C.140CCE 02 


VOLUME 


C.45456E 

05 

C.45319E 

05 

C.45181E 

05 

C.45044E 

05 

C.449C68 

05 

C.44769E 

05 

C.44632E 

05 

C.44494E 

05 

C.44357E 

05 

C.44219E 

05 

C.44082E 

05 

C.43945E 

05 

C.43807E 

05 

0.43670E 

05 

C.43532E 

05 

0.43395E 

05 

C.43258E 

05 

C.43120E 

05 

0.42983E 

05 

C.42845E 

05 

C.427C8E 

05 

C. 425716 

05 

0.42433E 

05 

C.42296E 

05 

0.42158E 

05 

C.42C21E 

05 

C.41884E 

05 

0.41746E 

05 

C.416C9E 

05 

C.41471E 

05 

C.41334E 

05 

C.41197E 

05 

C.41059E 

05 

C.4C922E 

05 

C.4C784E 

05 

0.40647E 

05 

C.40510E 

05 

C.4C372E 

05 

0.40235E 

05 

C.40097E 

05 

C.39960E 

05 

C.39823E 

05 

C.39685E 

05 

C.39548E 

05 

C • 394 10E 

05 

C.39273E 

05 

C.39135E 

05 

C.38998E 

05 

C.38861E 

05 

C.38723E 

C5 

0.38586E 

05 

C.38448E 

05 

C.3831 IE 

05 


DEL 


0.64790E 

00 

0.65868E 

00 

0.66943E 

00 

0.68014E 

00 

0.69082E 

00 

0.70147E 

00 

0.71207E 

00 

0.72265E 

00 

0.73318E 

00 

0.74368E 

00 

0.75415E 

06 

0.76458E 

00 

0.77497E 

00 

0.78533E 

00 

0.79565E 

00 

0.80594E 

00 

0.81619E 

00 

0.82641E 

00 

0.83659E 

00 

0.84674E 

00 

0.85685E 

00 

0. 86692 E 

00 

0. 87 696 E 

00 

0.88696E 

00 

0.89693E 

00 

0.90687E 

00 

C.91676E 

00 

0.92662E 

00 

C.93645E 

00 

0.94624E 

00 

0.95599E 

00 

C.96571E 

00 

C.97540E 

00 

0.98505E 

00 

0.99466E 

00 

0.10042E 

01 

0.10138E 

01 

0.10233E 

01 

0.10328E 

01 

0.10422E 

01 

0.10516E 

01 

0.10610E 

01 

0.10703E 

01 

0. 10796E 

01 

0 • 10888 E 

01 

0. 1098 IE 

01 

C.11072E 

01 

0.U164E 

01 

0.11255E 

01 

0. 11346E 

01 

0. 1 1436E 

01 

0.11526E 

01 

0.11616E 

01 
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